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Abstract 

Factorization of polynomials arises in numerous areas in symbolic com- 
putation. It is an important capability in many symbolic and algebraic 
computation. There are two type of factorization of polynomials. One is 
convention polynomial factorization, and the other approximate polyno- 
mial factorization. 

Conventional factorization algorithms use symbolic methods to get ex- 
act factors of a polynomial while approximate factorization algorithms use 
numerical methods to get approximate factors of a polynomial. Symbolic 
computation often confront intermediate expression swell problem, which 
lower the efficiency of factorization. The numerical computation is famous 
for its high efficiency, but it only gives approximate results. In this pa- 
per, we present an algorithm which use approximate method to get exact 
factors of a multivariate polynomial. Compared with other methods, this 
method has the numerical computation advantage of high efficiency for 
some class of polynomials with factors of lower degree. The experimental 
results show that the method is more efficient than factor in Maple 9.5 
for polynomials with more variables and higher degree. 

Key words: Factorization of multivariate polynomials, Interpolation 
methods, Numerical Computation, Decomposition of Affine Variety. 
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1 Introduction 



Polynomial factorization plays a significant role in many problems including the 
simplification, Grobner basis and solving polynomial equations etc. It has been 
studied for a long time and some high efficient algorithms have been proposed. 
There are two type of factorization of polynomials. One is convention polyno- 
mial factorization, and the other is approximate polynomial factorization. 

The modern conventional factorization methods follow Zassenhaus' approach 
[15] [16]. First, Multivariate polynomial factorization is reduced to bivariate fac- 
torization due to Bertini's theorem and hensel lifting[5][6]. Then one of the two 
remaining variables is specialized at random. The resulting univariate polyno- 
mial is factored and its factors are lifted up to a high enough precision. At last, 
the lifted factors are recombined to get the factors of the original polynomial. 

Approximate factorization is a natural extension of conventional polynomial 
factorization. It uses approximate methods to get approximate factorization 
of polynomial. The approximate factorization is not popular now, but there 
are some papers to discuss it. In 1985, Kaltonfen presented an algorithm for 
performing the absolute irreducible factorization, and suggested to perform his 
algorithm by floating-point numbers, then the factor obtained is an approxi- 
mate one. However, the concept of approximate factorization appeared first 
in a paper on control theory [TDj. The algorithm is as follows: At first express 
the two factors G and H of the polynomials F with unknown coefficients by 
fixing their terms, then determine the numerical coefficients so as to minimize 
\\F — GH\\. Huang et al. pursuit this approach, but the algorithm seems to 
be rarely successful, unless G or H is a polynomial of several terms. In 1991, 
Sasaki et al. proposed a modern algorithm)!!], which use power-series roots to 
find approximate factors. This algorithm is successful for polynomials of small 
degrees. Subsequently, Sasaki et al. presented another algorithm[12j which uti- 
lizes zero-sum relations. The zero-sum relations are quite effective for determin- 
ing approximate factors. However, computation based on zero-sum relations 
is practically very time-consuming. In [13] . Sasaki, T. presented an effective 
method to get as many zero-sum relations as possible by matrix operations so 
that approximate factorization algorithm is improved. In [3], Corless et al pro- 
posed an algorithm for factoring bivariate approximate polynomial based on 
the idea of decomposition of affine variety. However, it is not easy to generalize 
the algorithm to factor multivariate approximate polynomials. Recently, Zhang 
et al|17j proposed an algorithm for obtaining exact rational number from its 
approximate floating number. In this paper, basing on the algorithm in [17], 
we present an algorithm which use approximate method to get exact factors 
of a polynomial. It can be regard as a generalization of Corless' algorithm in 
multivariate polynomial and exact polynomial case. 

The remainder of the paper is organized as follows. Section 2 presents a for- 
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inula by which a polynomial is constructed from sampled points on its variety. 
A condition is given for the formula to determine only one polynomial up to 
a nonzero constant multiples, and the error estimation is discussed. Section 3 
gives a review of a modified continued fraction method, by which an exact ratio- 
nal number can be obtained from its approximation. Section 4 first discusses the 
error control, and then proposes a factorization algorithm for multivariate poly- 
nomial over rational number field. Section 5 gives some experimental results. 
The final section makes conclusion. 



2 Interpolation method 

Polynomial interpolation is a classical numerical method. It is studied very well 
for univariate polynomials. In general, there are four types of polynomial in- 
terpolation method: Lagrange Interpolation, Neville's Interpolation, Newton's 
Interpolation and Hermite Interpolation. Lagrange interpolation formula can 
get the interpolation polynomial at once for a given set of distinct interpolation 
points and corresponding values (xi, /,), i = 1, ■ • • , n+1. It is very useful in some 
situations in which many interpolation problems are to be solved for the same 
set of interpolation points Xi, i = 1, • • • , n+1, but different sets of function values 
fi, i = 1, • • • , n + 1. Unlike Lagrange interpolation method which solve the inter- 
polation problem all at once, Neville's interpolation method solve the problem 
for smaller sets of interpolation points first and then update these solutions to 
obtain the solution to the full interpolation problem. It aims at determining the 
value of the interpolating polynomial at some point. It is less suited for deter- 
mining the interpolating polynomial. If the interpolating polynomial is needed, 
Newton's Interpolation formula is preferred. Just like Neville's interpolation 
method, it first get interpolating polynomial for smaller sets of interpolation 
points and then update the polynomial for a larger sets of interpolation points, 
step by step, and finally, the interpolating polynomial is obtained for the set 
of the full interpolation points. If the interpolating problem prescribes at each 
interpolation point Xi,i = 1, ■ ■ ■ ,n+l not only the value but also the derivatives 
of desired polynomial, then the Hermite formula is preferred. 

For univariate polynomial interpolation, n+1 distinct interpolation points 
and their corresponding values determine only one polynomial with degree less 
than or equal to n. However, the interpolation points for multivariate poly- 
nomial interpolation can not be chosen arbitrarily. They need to satisfy some 
conditions. So we need a definition as follows: 

Definition 1 Let be a set of n-dimension points and P a polynomial space. 
We call Proper interpolation points of P if for any f defined on 0, there 
is a unique polynomial p E P matching f at 0. 
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In definition [TJ a polynomial p matching / at means p\Q — f\Q. In 
general, we can determine interpolation polynomial space such as: if knowing 
the total degree d of /, we choose P = {p\deg(p) < d ,p £ K\x\,x 2 , ■ ■ ■ ,x n ]}; 
if knowing the degree dj of / in n(i = 1, • ■ • , n), we choose P = {p\ /\t=i,i,—,n 
{deg Xi (p) < di),p € K[x\, X2, ■ ■ ■ , x n ]}. Once the interpolation polynomial space 
is determined, the proper interpolation points O of P can be set by interpolation 
methods [9] pQ. 

In this paper, we need to construct a polynomial from some points of its 
variety. Values of the polynomial at interpolation points are all zero. So, we 
introduce an interpolation formula for this case. 

Let f(x\, X2, • • • , x n ) be a polynomial to be interpolated. It is represented as 
follows: 

f(x u x 2 , ■ ■ ■ , x n ) = Cl X a > + c 2 X a > + ■■■ + c m X am , (1) 

where X ai the distinct monomials, and Cj are the cor- 

responding coefficients. 

Let Pi,P2, ■ ■ ■ iVm— 1 be points on variety of f(xi,X2,---,x n ), where = 
(Pn,Pi2, ■ • ■ ,Pin) ior i = ,m - 1. P" J = p^p^' 1 ■ ■ -pfn' 3 denote the 
value of the monomial X aj at Pi. An interpolation formula is as follows: 
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Next, we need to know what condition the interpolation points pi,P2, • • • ,Pm~i 
should satisfy so as to ensure formula @ to determine a unique polynomial and 
it is f(xi, X2, ■ ■ ■) up to a nonzero constant multiplies. 

Let X* denote the minor of X ai resulting from the deletion of row 1 and 
column i in formula @. We have a theorem as follows. 

Theorem 1 Let f(x\,X2, ■ ■ ■ ,x n ) be a nonzero polynomial and it is expressed 
as in equation (OJ). If the m—1 zeroes of f(xi, x 2 , • • • , x n ) satisfy a condition that 
G(xi, • ■ • , x n ) 7^ in formula^ then formula ^determines a unique polynomial 
and it is polynomial f(xi, X2, ■ ■ ■ , x n ) up to a nonzero constant multiples. 

Proof: Due to f(x\, X2, ■ • ■ ,x n ) ^ 0, the coefficients of f(x\, X2, ■ • ■ ,x n ) are 
not all zero. Assume that ^ 0, Cj 2 ^ 0, • ■ • , Cj s 7^ 0, and their corresponding 
monomials are X ai i , X ai ? , • • • , X ai ° . Let X* , X* 2 , • • • , X* denote the minors 
of A tti i , A Qi 2 , • • • , X ai * in formula respectively. 

First, we show that if one of X* ly X* 2 , • • • , X* s is nonzero, then f{xi,x 2 ,--- 1 x n ) = 
cG(xi, X2, ■ ■ • , x n ), where c is a nonzero constant. Without loss of generality, let 
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us assume that c k ^ and X£ ^ 0. So we have 
G(x L ,x 2 , ■ ■ ■ ,x n ) 
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For i = 1, 2, ■ 
yields: 



■ ,k — 1, fe + 1, • • ■ , m, addition of Cj times column i to column A; 



G(x 1 ,x 2 , ■ ■ ■ ,x n ) 
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Due to Ck ^ and Xj! 7^ 0, it follows that f(xi, X2, • • ■ , x n ) = cG(x\,X2, ■ ■ ■ ,x n ), 
where c = Ck is nonzero constant. 

Second, we assert that one of X* l7 X* 2 , ■ ■ ■ , X* d must be nonzero. We prove 
it by contradiction. Let us assume that X* x = 0, X* 2 = 0, • • • , X* d = 0. 

Under the assumption of the theorem, G(x\, X2, x n ) = Y^iLi biX ai ^ 
and f(x\, X2, • • • , x n ) ^ 0. So, not all of their coefficients are zero. Assume 
that b hl ^ 0, b h2 ^ 0, • • • , b hz ^ and c ix ^ 0, c i2 ^ 0, • • • , c is ^ 0. Since 
X* x = 0,X* 2 = 0,---,X* s = 0, it holds that i k ^ hj for k = 1, • • • , s and 
j = 1, • • • , z. Hence, we have that X ai k ^ X ah i for k = 1, • • • , s and j = 1, • • • , z. 
Let 

H(xi,X 2 , ■ ■ ■ ,X n ) = f(xi,X 2 ,---,Xn) + G(xi,X2,---,X n ) 

= Ql X a " + • • • + c ie X ai ° + b hl X a ^ +■■■ + b h X a ^ 
+ 
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Because of ^ and X^ ^ 0, it has been shown above that 

G(xi, x 2 , ■ ■ ■ , x n ) = cH(xi,x 2 , • • • , x n ), 

where c is nonzero constant. Hence we deduce that 

G{xi,x 2 , ■ ■ ■ ,x n ) ~ cH{xi,x 2 , ■ ■ ■ ,x n ) = 

. However it is impossible because the term cc^X ^ ^ is not monomial of 
G{x\ , ■ ■ ■ , x n ) . Therefore we show that one of X* x , X* 2 , • • • , X* s must be nonzero. 
The proof of the theorem is finished. 

The above theorem shows that if formula @ gives a nonzero polynomial, and 
the first row contains all monomials of an interpolating polynomial, then it is 
the interpolating polynomial up to nonzero constant multiples. 

However, due to floating-point computation, we only get approximate zeroes 
of /. Accordingly, we only obtain approximate factors. In the remaining of this 
section, we study how the error resulting from floating computation affects that 
of factors. For simplicity, let us give a definition: 



Definition 2 Let X™ 3 = xTTxfT ' ' ' X TT > where x^i, . . . , Xi >n are complex 
number and a^i, • • • , ctj, n are nonnegative integer. A generalized Vandemonder 
determinant is defined as follows: 
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We have an estimation of generalized vandemonder determinant as follows: 



Theorem 2 Let M = max^jX^} and B = m&x itjik {\\XT - Xf h \\, \XT 
XT II}- Then for m > 2 it holds that 

\V m \<m\M m - l B 



Proof: We prove it by inductive method. When m = 2, the generalized Van- 
demonder determinant is 
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Assume that \V m \ < m\M m 1 B for m — k. Let us show that it holds for 
m = k + 1. We expand Vk+i by minors as follows 

\v k+ i\ = \J2(- 1 ) i+jx i Jv k j \ < E \ x ?\* K j \ 
3=1 3=1 

According to our assumption that V fc iJ < k\M k ~ 1 B, we have 

fc+i fe+i 

|Vfc+i| < \ X t 3 1 * K' 3 \ < M * klMk ~ lB = ( k + l) } -M k B 
3=1 3=1 

The proof is finished. 



Theorem 3 Let M = max^Xf 3 } and B = max iJjfe {||Xf' - (Xj** - 

A^ l ||} and e = max^j |ai|. ^4 determinant is as follows: 
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Then we have an estimate that \V m \ < M m 2 mlBe for m > 3. 



Proof. Expanding V m by column (oi, • • • , a m ) T and then using theorem^ we 
can get the proof. 

And now, we study the difference between two generalized Vandemonder 
determinants. 
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Proof: 
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From theorem [3j it holds that \Vm^ 
The proof is finished. 
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3 Continued fraction method 



As we said above, our method is to use approximate method to get exact factors 
of a multivariate polynomial over rational number field. So we need to recover 
the exact coefficients of a polynomial from its approximate coefficients. In this 
section, we introduce a continued fraction method to recover exact rational num- 
ber from its approximation. As we know, a continued fraction representation of 
a real number is one of the forms: 

fl o+ , 1 i , (4) 

°i + - i i 

where ao is a integer and ai, a 2 , 03, • • • are positive integers. One can abbreviate 
the above continued fraction as 

[a ;ai,a 2 , ■ ■ ■} 



In order to recover exact rational number, we introduce a control error into 
the conventional continued fraction method. The continued fraction method is 
modified as follows. 
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Algorithm 1 Continued fraction method 

Input: a nonnegative floating-point number a and Z\ > 0; 

Output: a rational number b 

Step 1: i := 1 and x\ := a; 

Step 2: Getting integral part of Xi and assigning it to di, assigning 
its remains to bi . If bi < e±, then goto Step 

Step 3: i:=i + I; 

Step 4- x i '■— an d goto Step\3[ 

Step 5: Computing expression and assigning it to b. 
Step 6: return b. 

In [17], we discussed how to get error control E\ > 0. The theorem is as 
follows: 

Theorem 5 Let no/n\ be a reduced rational number and r its approximation. 
Assume that nQ,n\ are positive integers and L > max{ni, 2}. K is a positive in- 
teger. The continued fraction representations of n^/ni and r are [do, a%, ■ ■ ■ , a/\r] 
and [bo, &i, • ■ • , 6m] respectively. If \d\ = \r — n /ni\ < 1/((2K + 2)L{L — 1)), 
then one of the two statement holds 

• ai = bi for i = 0, • • • , N , and b^+i > K; 

• ai = bi for i = 0, ■ • • , N — 1, and b^ = — 1. &jv+i = 1, > K- 

From theorem [5j getting exact non-negative number 712/ni from its approx- 
imation ro is summarized as follows: 



Algorithm 2 Obtaining Exact Number 



Step 1: estimating an upper bound of the denominator 0/712/^1, De- 
noted by L; 

Step 2: computing 





(2L + 2)L(L - 1)) 



Step 3: obtaining r Q by approximate method such that \r — n 2 /ni\ < 

(3; 

Step 4 : taking e\ = 1/L in algorithm^ and calling algorithm\]\to get 
b. So n^/rii = b. 
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4 Factoring Multivariate Polynomials by Approx- 
imate Method 



In this paper, we only discuss factorization of a multivariate polynomial over 
rational number field. So its coefficients are all rational numbers. In order 
to get factors of a multivariate polynomial over rational number field, we first 
compute its factors over complex field. These factors are complex coefficient 
polynomials. Products of some of them must be real polynomials. We get 
the products which are approximate rational coefficient factors of the original 
polynomial. Finally, transforming these real products into rational coefficient 
polynomials yields factors of the original polynomial over rational number field. 

A set 

V(f) = {(oi, • • • , O G C" : /(oi, ■ • • , a n ) = 0} 

is called affine variety of f(xi, ■ ■ ■ , x n ). An affine variety V C C n is irre- 
ducible if whenever V is written in the form V = V\ U V2, where V\ and V2 are 
affine varieties, then either V\ = V or V% = V. 

Let f(%i, X2, • ■ • , x n ) be a square free polynomial over complex number field 
C and / = /1/2 ■ • • f m , where /j is distinct irreducible polynomials. Then < / > 
is a radical ideal. It holds as follows 

V(f)=V(f 1 )UV(f 2 )U---UV(f m ), (5) 

where V(fi) are irreducible affine varieties. 

From equation ([5]), if we get a point on variety of /, it must be either on 
one of V(fi) or on the intersection of them. When the point is not singular 
point, it must be on one of V(fi) and not on the intersection of two varieties. 
Theorem [T] shows that if getting enough points in some variety of V(fi) that 
satisfy the condition of theorem [TJ we can recover the polynomial by formula 
([2]). Therefore, the procedure of factorization is as follows: First get a initial 
nonsingular point on one variety of V(fi). And then obtain enough sampled 
points on the same variety. Third, use formula {2j to get a factor and finally, 
obtain a rational factor. 

However, due to approximate computation, we first discuss error control, and 
then study factorization. 

4.1 Error control 

Let f(xx, ■ ■ ■ , x n ) be a polynomial to be factored over rational number field. 
According to algorithm [31 the first thing we need to do is to determine an upper 
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bound of absolute values of coefficient denominators of factors of f(x%, X2, • • • , x n ). 
The following theorem is very helpful. 

Theorem 6 Let g(xi, ■ ■ • , x n ) be a monic polynomial over rational number 
field. Its factorization over rational number field is g(x\, • • • , x n ) = Yii=i hi{x\, ■ ■ • 
where all hi{x\, • • • , x n ) are monic polynomials. Assume that N > is the least 
common multiple of denominators of coefficients of g(x\, • ■ • , x n ). Then N is an 
upper bound of absolute values of denominators of coefficients of hi(x%, • ■ • , x n ) 
for i = 1 , • • • , m . 

Proof: It is clear that Ng(x%, ■ ■ ■ , xjv) is a primitive integral coefficient poly- 
nomial. Let Ni be the least common multiple of denominators of coefficients of 
hi(xi, ■ ■ ■ ,x n ). Hence, hi(xi, ■ ■ ■ , x n ) = ■j^h i (xi, ■ ■ ■ ,x„),where hi{x\, ■ ■ ■ , x n ) 
is a primitive integral coefficient polynomial. From Gauss'lemma, YIaLi h%{x\, • • ■ , 
is a primitive polynomial over Z[x\, • • ■ , x n ]. On the other hand, we have that 

m m 

Ng(xi,---,x n ) = NY\_hi(xx, - ■ ■ ,x n ) = =fm — JJ hi(x%, ■ ■ ■ ,x n ) 

Since Ng{x\, • • ■ , x n ) and YiiLi hi(xi, • • ■ , x n ) are primitive polynomials, it holds 
that N/^XiL\ N i) = ±l ■ Therefore, N = ± J]™ i N i- Th e proof is finished. 

Theorem [6] shows that the positive least common multiple of denominators 
of coefficients of a monic polynomial is also an upper bound of absolute values 
of denominators of coefficients of its monic factors. 

According to algorithm[2]and theoremlU we calculate control error as follows: 



Algorithm 3 Calculating control error 

Step 1: calculating upper bound of absolute values of the coefficient 
denominators of exact factors, denoted by L. From theorem^ we 
can take the positive least common multiple of the coefficient denom- 
inators of a monic polynomial to be factored. 

Step 2: taking K = L + 1 and E\ = 1/K in algorithm^ 
Step 3: computing (3 = (2K+2)l(l-i)) ■ 

Step 4 : computing the control errors in theorem^such that mm\M m ~ 2 Be < 
(3. 

Therefore, in order to factor polynomial f{x±, • • ■ ,x n ), we first call algorithm [3] 
to compute control error e and control error e\ in algorithm [T] 
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4.2 Initial point 



Let f{x\, x 2 , ■ ■ ■ , x n ) be a square-free polynomial to be factored over rational 
number field. Of course, it is a square-free polynomial over complex number 
field. Choosing n — 1 floating-point numbers Xi t o, £2,0, • • • , x ra -i,o at random 
and numerically solve f{xi$, X2,o, ■ ■ ■ , x n -i,o, x n ) = for variable x n within 
control error e. Denote this solution by x n ,o- So the zero (£1,0, £2,0, ■ • • > Xn,o) 
of f(xi, ■ • ■ , x n ) must be in one of V(fi) or the intersection of them. If 

V/(xi,o, x 2 ,o, x n ,x) (0,0, • • • , 0), 

then it must be in one of V(fi) and not in the intersection of two of them. 
In order to get a neighborhood Uq of (xi.o, £2,0, • ■ • , x n fi) which is only on one 
variety V(/j), we require that 



- + 



(6) 



9/(^1, 0,"-,^m,o) / Q 

3a;„_i ' 



We always assume that anx £: ,Xn) is not zero polynomial. This is because 
if there exists fc(0 < k < n) such that 9f{xi Q^ k "' Xn) = 0, then 

f(xi, • • • , x k -i,x k ,Xk+i,- ■ ■ , x n ) = g(xi, • • • , x k -i,Xk+i,- ■ ■ ,x n ). 

So, we can consider polynomial g(x\, • • ■ , Xk-i, JEfe+i, • • • , x n ). Hence with prob- 
ability 1, we can get a point (xi.o, X2,o> • • * ■> x n ,o) that satisfies equation©. 



4.3 Sampled points 



Let X = (xi,X2, ■ ■ ■ ,x n ) and P = (xx,o,x 2 ,o, ■ ■ ■ ,x n>0 ). Since 'q^' x ^ 
for i = 1, • • • ,n — 1, we can calculate h± > 0, h 2 > 0, • • • , /i n -i > from Mean 
value theorem such that 

df df df df 

{sign{ — {P )), ■ ■ • , stgn(— — -(P ))) = ( si 5 n (^"( X ))' - ' ' ' Slgn (fa — ~ W)) 

for X 6 f(Po), where signQ is sign function and let U{Pq) = [[2^1,0 — hi,x%fi + 
h{\, ■ ■ ■ , [x n _i,o — h n -x, j; n -i,o + ^n-i]J- 

cZf is denoted by the degree of polynomial /(xi, ■ ■ ■ , x n ) with respect to Xi for 
i = 1, • • • , n — 1. Choosing di + 1 distinct points in interval [x^o — ^i, x^o + hi] for 
i = l,2,---,n — 1. So we get n— 1 dimension vectors denoted by And then 

substituting each vector € $> n -i into f{xi, x 2 , ■ • • , x„) and getting univariate 
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polynomial g(x n ), Solving g(x n ) = within control error e in algorithm [3] and 
choosing the solution x n ^ which is the closest to x n ,o yields a n-dimension vector. 
Hence we get Yli=i (^» + 1) ^-dimension vectors denoted by 4>„. 

4.4 Getting a factor by sampled points 

In this subsection, we discuss how to get exact factors of a polynomial. This 
procedure runs as follows: First a candidate set monomials in the support of 
the factor is selected. If we know the pattern of monomials of the factor, we 
choose a restricted set of monomials. If we do not know the pattern, then 
for m = 1, • • • , d — 1, where d is total degree of f{x\,X2, • • • , x n ), we use the 
complete set of monomials of degree less than or equal to to, denoted by M m . 
Next, for every set of monomials M m (m — 1, 2, • • • , d — 1), some of the above 
sampled points are selected. If the sampled points is less than we want, we can 
refine U (Pq) and get enough sampled points. The selected sampled points must 
satisfy the condition of theorem [TJ From the selected sampled points, formula 
^ is used to obtain a polynomial g{x\, X2, ■ ■ • , x n ) which monomial set is M m . 
and then transform g(xi, x%, ■ ■ ■ , x n ) into a monic polynomial, still denoted by 
g(x\,X2, ■ ■ ■ ,x n ). Third, we deal with g(x±,X2, ■ • ■ ,x n ) in two cases: 
Case 1: g(xi,X2, • • ■ ,x n ) is not a real polynomial but a complex polynomial. 
Let the selected sampled points be po,Pi, • • • , from which g(x\, • • • , x n ) is con- 
structed. So, their complex conjugate points po,J>i, ■ • •, are on the other variety 
V{fj). We construct a monic polynomial g{x\, x%, • • • , x n ) from po^Pi, 1 • •• So 
gg is a monic real polynomial. We deal with gg just as in case 2. 
Case 2: g{x\, ■ ■ ■ ,x n ) is a monic real polynomial. We use the polynomial di- 
vision proposed in [3] to get h{x\,X2, ■■■ ,x n ) which minimizes ||/ — gh\\. Let 
r = \\f — gh\\. If r is large, we should add more sampled points and extend 
the monomial set M m to M m +i, and then use formula ((2|) to get a poly- 
nomial with higher degree. If r is very small, then we use algorithm [5] to 
transform g(x±, ■ ■ ■ ,x n ),h(xi, ■ • ■ ,x n ) to a rational polynomial g'(xi, ■ • ■ , x n ) 
and h'(xi, ■ ■ ■ , x n ) respectively. Compute r' = \\f — g'h'\\. If r' = then 
g'{xi, • • • , x n ) is a rational factor of f{x\, • • • , x n ). We continue to factor poly- 
nomial h'(xi, • • • , x n ) on the other variety V(fj) . If r' ^ 0, then from theorem 
g(xi,X2, ■ ■ ■ ,x n ) is not an approximate rational factor but an approximate 
real factor of polynomial f(xi, x%, ■ ■ ■ , x n ). We continue to factor polynomial 
h(xi, • • • , x n ) on the other V(fj). Finally, we get all factors of f(x\, • ■ • , x n ). 
Let f(xi, ■ ■ ■ ,x n ) = g\g2 ■ ■ ■ gugk+i ■ ■ ■, and $ = {31,52, • ■ • ,gk} are not rational 
factors, the others are rational factors. We compute products of two distinct 
polynomials of $: = gi * gj(i ^ j,and i,j < k). Check if every gij is an 
approximate rational factor, Whenever it is so, keep g^ as a rational factor 
and remove gi,gj from <E>; After finishing to deal with products of two distinct 
polynomials in $, we compute products of three distinct polynomials of <f>: 
gijm — gigjgm {i 7^ h3 7^ m,m ^ i), and check if every g ijm is an approximate 
rational factor. Whenever gij m is so, keep it as a rational factor and remove 
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9ii9j,9m from $; and so on, until all rational factors of f(xi,x 2 , ■ ■ ■ ,x n ) are 
obtained. 



5 Experimental results 

The following five polynomials are randomly generated by randpoly command 
in Maple. Our algorithm is implemented in Maple. Compared our algorithm 
with factor command in Maple in the platform of Maple 9.5 and PHI LOG, 
256M RAM, The running time of the five examples are as follows: 

Example 1. This polynomial is with four variables, 46376 terms, and of 
degree 30. The 15 factors are of degree 2. The factor command in Maple costs 
430.108 seconds, our algorithm costs 278.375 seconds. 

Example 2. The polynomial is with four variables, 52360 terms, and of 
degree 31. The 21 factors are of degree 1 or 2. The factor command in Maple 
costs 667.438 seconds to factor the polynomial, our algorithm 307.922 seconds. 

Example 3. The polynomial is with four variables, 52360 terms, and of 
degree 32. The 32 factors are of degree 1. The factor command in Maple 
computed 7200 seconds and gave no result, our algorithm 250.436 seconds. 

Example 4. The polynomial is with three variables, 17296 terms, and of 
degree 45. The 15 factors are of degree 3. The factor command in Maple use 
261.468 seconds to factor the polynomial, our algorithm 173.265 seconds. 

Example 5. The polynomial is with three variables, 37820 terms, and of 
degree 59. The 25 factors are of degree 2 or 3. The factor command in Maple 
used 1683.655 seconds to factor the polynomial, our algorithm used 710.045 
seconds. 



6 Conclusion 

• Our algorithm first need a initial point. We just get n — 1 floating-point 
number at random, and then substitute them into a polynomial to be fac- 
tored and obtain a univariate polynomial. Solving the univariate polyno- 
mial for the last variable x n yielding d n solutions. In the above discussion, 
it seems only to take one solution and throw away the other solution. In 
fact, we should keep these d n — 1 points as initial points on the other vari- 
eties. For the same reason, we should keep other d n — 1 solution whenever 
get a sampled point in neighborhood Uq. Once we get sampled points in 
neighborhood Uq, we obtain d n — 1 sets of sampled points on the other 
d n — 1 varieties of the original polynomial. 
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• Our algorithm is to get exact factorization of polynomials by approximate 
method, so its efficiency is higher than symbolic factorization method 
when a polynomial is with more variables and higher degree. However 
its efficiency is lower than symbolic factorization when a polynomial is 
with less variables and lower degree. In order to take advantage of both 
numerical and symbolic factorization algorithm, we can factor a polyno- 
mial with more variables and higher degree as follows: First, we use our 
algorithm to get some factors of the polynomial and remove these factors 
from the polynomial. When the remaining polynomial is with less vari- 
ables and lower degree, then, we use symbolic factorization algorithm to 
get the remaining factors. In fact, Our algorithm is implemented by this 
idea. 

• Noting that interpolating formula ^ is lower efficient when the factors 
of the original are with more variables and higher degree. So we should 
improve interpolating formula Q further. 
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